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SUMMARY* 


The first part of this paper pertains to the estimation of subsonic 
rotary stability derivatives of wings. The unsteady potential flow 
problem is solved by a superposition of steady flow solutions. Numerical 
results for the damping coefficients of triangular wings are presented 
as functions of aspect ratio and Mach number, and are compared with 
experimental results over the Mach number range 0 to 1* 

In the second part, experimental results are used to point out a 
close correlation between the nonlinear variations with angle of attack 
of the static pitching -moment curve slope and the damping -in-pitch 
coefficient. The underlying basis for the correlation is found as a 
result of an analysis in which the indicial function concept and the 
principle of superposition are adapted to apply to the nonlinear problem. 
The form of the result suggests a method of estimating nonlinear damping 
coefficients from results of static wind-tunnel measurements. 
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ESTIMATION OF ROTAKI STABILITY DERIVATIVES 


AT SUBSONIC AND TRANSONIC SPEEDS 

By Murray Tobak* and Henry C. Lessing* 

National Aeronautics and Space Administration 
Ames Research Center 
Moffett Field, Calif. 


1. INTRODUCTION 

It is now generally recognized that modem aircraft, particularly 
tailless aircraft, can experience a significant loss of damping in the 
short-period pitching mode at transonic speeds. This loss has been 
traced to a reduction or even a change in sign of the damping -in -pitch 
parameter + C^. 0 n the supersonic side of the transonic speed 

range a fairly complete understanding of the mechanism underlying the 
reduction of pitch damping has been made possible by the rapid develop- 
ment of theory in this range (cf. (l) and attendant bibliography). 

It is known, however, that the phenomena actually have their origins in 
the subsonic speed range; unfortunately, the great difficulty of the 
subsonic theoretical problem for finite -span wings has prevented a 
detailed study of these origins. Moreover, in the transonic speed range 
itself, and in fact, at all speeds, when flight conditions depart from 
those implicit in the range of applicability of a linearized theory, it 
is known that the aerodynamic forces and moments can become highly non- 
linear functions of angle of attack. Here again, the great difficulty 
of the theoretical problem has prevented our gaining an understanding 
of pitch-damping behavior under such circumstances. While it may not be 
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feasible to obtain explicit theoretical solutions in these cases,, the 
possibility remains to formulate the problems so they relate closely to 
analogous problems in steady flow. The aerodynamics of steady flows 
being of fundamental importance in aircraft design, there exists a far 
greater fund of knowledge, both theoretical and experimental, relating 
to this field than exists for unsteady flows. Hence, if it is possible 
to relate an unsteady flow problem to an analogous problem in steady 
flow, one has the hope of bringing this fund of knowledge to bear. 

This possibility will be examined for the two problems mentioned above; 
first, in the estimation by linear theory of pitch-damping coefficients 
of wings flying at subsonic speeds; second, in the estimation of pitch- 
damping coefficients when the aerodynamic forces and moments are non- 
linear functions of angle of attack. 

2. ESTIMATION OF Cj. AND FOR WINGS AT SUBSONIC SPEEDS 

Consider the first problem, the estimation of pitch damping of 
wings flying at subsonic speeds. As is well known, to a first order in 
frequency the pitch-damping coefficient is composed of two parameters; 

Cniq, the pitching-moment coefficient proportional to constant pitching 
velocity, and C^, the pit ching -moment coefficient proportional to 
constant vertical acceleration.- Now the analysis of C^ already 
corresponds to that of a steady flow, namely the flow over a wing cambered 
and twisted to have a downwash distribution linearly dependent on chord- 
wise distance, x. Hence, all attention can be focused on the remaining 
parameter, C “d* 
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2.1 Solution to Potential Equation 

To reduce the calculation for to a steady flow problem, it is 

possible to adapt a simple device previously used extensively in calcu- 
lating Cm£ at supersonic speed. The device was first introduced by 
Ribner and Malvestuto (2), who in turn credit it to C. S. Gardner. 

A solution is sought to the unsteady wave equation 


3 cp. 


•XX + V + 


Pzz ~ a^s ^xt “ a^2 ^tt 0 


( 1 ) 


subject to the boundary conditions in the z = 0 plane 
w(x,y, o) = cp z (x,y,o) = -aVt , on the wing 

^ f (*x + k 9 t) z=o = 0 ' off the ^ 


( 2 ) 


Gardner has shown that a solution to Eq. (l) suitable for supersonic 
speed can be built of a combination of steady flow potentials; namely the 
potential due to steady unit pitching velocity i(r, and the potential due 
to steady unit angle of attack X, in the combination shown in Eq. (3) 

cp(x,y,z,t) = - £2 ^ + a (t + X a=1 (x,y,z) 

+ <i> c (x,y,z) . ( 3 ) 

The adaptation of Eq. (3) to make it applicable to subsonic speeds consists 
only in adding a third steady flow potential, <J>c • The necessity of this 
addition will be evident from inspection of the loading coefficient 
expression. 


2.2 Loading Coefficient 

The loading on the wing is obtained from Bernoulli's equation 


4p _ 4 
V 


<P x (x,y,o,o) 


+ i cp t (x,y, 



( 4 ) 
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and this leads to 

+ (0 C (5) 
evaluated at t = 0. Now if the flight speed is supersonic and the wing 
has supersonic trailing edges, the wake is of no concern in calculating 
loading on the wing since it cannot affect conditions on the wing . In 
these circumstances, the first three terms in Eq. ( 5 ) constitute the 
solution for loading on the wing. For subsonic speeds, however, the 
Kutta condition must be invoked; this entails that loading be zero at the 
wing trailing edge, and of course, also in the wake. This condition is 
fulfilled by the first two terms in Eq. ( 5 ) but not by X itself. 

Writing X as 



one sees that X is a function of y in the wake, retaining at each 
station y the value it has at the wing trailing edge. This value is 
directly proportional to the span loading due to angle of attack, which 
may be presumed to be known. Hence, the potential is added so that 

its loading contribution will cancel a spurious but known loading in the 
wake contributed by X. In addition, since the boundary condition for 
downwash on the wing is already satisfied by i|r and X alone, a second 
condition to be fulfilled by $ c is that 3$ c /dz be zero at the wing 
surface. These two conditions suffice to" determine the additional loading. 

2-3 Solution for Triangular Wing 

To avoid excessive mathematical development, it is merely stated here 
that after some manipulation, the problem for the correction loading can 
be reduced to finding the loading corresponding to a known downwash 


Ap _ acp 

n v 


m 2 f Ap\ • | M 2 x Ap\ 
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dependent only on spamri.se distance, y. Then by use of the reverse flow 
theorem, the lift and moment contributions are found relatively simply 
by spanwise integrations of this dovnwash distribution, multiplied by the 
reverse-flow span loadings due to angle of attack and pitching velocity. 
The form of the result as it applies to the triangular wing is given 
in Eq. (7) . 

= Cl^ + -y" 1 g(u)L(u) R ^ a=1 du 


■ So + / 1|2a£ (5 ’’ l) F,a= 1 d ’> +/ 1 g(u)«u) 


> (7) 


where 


g(u) = ^ r Ap 


(I >4^,0=! 


dq 


s/(l -|) 2 + P 2 m 2 (u-r|) ; 


7.| <3co 

(t ' u) S,a=l« 

L(u) Rjq=1 = jj^| ^ ^ ’ U ^R, q=i d £ > 3x13 at a P ex 


and the subscript o on the stability derivatives is meant to indicate 
that the terms are to be evaluated for an axis located at the wing apex. 
The subscripts F and R distinguish between loadings in forward and 
reverse flow, respectively. 


2.4 Approximate Steady-State Loadings 

In the form given, Eq. ( 7 ) is exact. It will be noted that the 
calculations involve the steady-state loadings due to angle of attack 
and pitching velocity, and the corresponding span loadings in reverse 
flow. Hence, the stated objective has been achieved. To complete the 
calculations, it would of course be most advisable to use experimental 
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data for these quantities, or secondly, exact theoretical solutions. 


Since probably, neither -will be available in sufficient detail, approxi- " 
mat ions must be introduced at this stage. The calculation : is a sensitive 
one, especially for pitching moment, in that it involves differences 
of numbers of approximately equal size. Therefore, in introducing 
approximate loadings it is advisable to use a set of loadings that is 
at least internally consistent, that is, a set that satisfies the reverse- 
flow theorem. In calculations for the triangular wing, it is possible 
to achieve this aim by a generalization of a technique introduced by 
Lomax and Sluder (3); namely, it is assumed that the loadings are those 
given by slender -wing theory, multiplied by a chordwise correction factor. 
The correction factor alters the slender-wing loadings both to account 
for compressibility and to comply with the Kutta condition. The loadings 
have the form given in Eq. (8) . 


wp = w F (x) 
w R = v R (x) 



It will be noted that by virtue of the form of Eqs . (8), the integrand 
of the reverse-flow relation will form a perfect differential. The 
integration in x over the chord is then identically zero since f is , 
zero at the trailing edge and the square root is zero at the leading 
edge. Hence satisfaction of the reverse-flow theorem is insured in a 
very simple way. The • chordwise correction factor f(x/co) may be obtained 
from results presented in (3)* •. 
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Integrating the first of Eqs. (8) to obtain the steady-state 
stability derivatives in forward flow yields 

Cj,, . = -jcA /' 1 | 2 f(|)d.| 

J ° x > (9) 

= | ^/^ 2 f(l)d| 

' °\,’ -i”"X 1?3f(e)d| 

Note that by virtue of Eqs. (8), md'^ are related simply by the 

factor 3/2. It has also been found that the center of pressure of the 
lift due to pitching about the apex is given very closely by 9/8 times the 
center of pressure of the lift due to angle of attack}* that is 

°\, = ~ 2 

C mq^ _ 9 J (l°) 

C 'Lq 0 '8 Cla 

Inserting Eqs. (8), (9) and (10) in Eqs. (7) leads to the following 
expressions for Cj^, and 

= | + °L (g) 

( 11 ) 

' C n& 0 = " ^ ^a 2C Ia + 


where 


^l(s) = - / g(u)L(u) R ^du 
J o 

Cmjg) =J^ 1 g(^)L(^) R ,q= 1 , ap ex du 

Curves of the parameters Cx w /( rtA./2 ) , Xq,, and the correction terms 


*This result is also obtained using Lawrence’s method (4) of computing 
the steady-state loadings 
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0 2 C L (g)/(iiA/2), P 2 C mo (g)/(j tA/2) are presented on Figs. 1 and 2. 

From these results, all the stability derivatives can be obtained for 
triangular wings within the reduced aspect ratio range 0 < PA < 4. 


2.5 Empirical Corrections 

Observe that by virtue of Eqs. (lO) the expressions for the stability 
derivatives have been put almost completely in terms of Cj^ and Xq,. It 
is to be expected that an improvement in accuracy and, for example, at 
least a partial account of wing-body interference effects could be 
realized by replacing Cj^ and x a with values determined from experi- 
ments. The given results can be corrected in the following way. Let 




°S h 


^exp 



j 


(12) 


Then 


So ■ I 

So " ' § 


S- “ 


/ (13) 


2 ^h 0 ^ + ^®)thj . 

_ k v Sh a °S h + “S>(e)th] 

These are the results which will be compared with results of experiments. 


= [IV 


2.6 Comparison with Experiments 

The theory as computed and corrected from Eqs. (13) is compared with 
experimental data for triangular wings of aspect ratio 1.45, 2, and 4 on 
Figs. 3 to 5* Experimental results for the A = 1.45 wing were, obtained. 
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at the Swedish Aeronautical Research Institute ( 5 ); those for the A = 2 

t 

and 4 wings, at the Ames Research Center. (6), (j). On each figure, the 
damping coefficient as defined here is shown multiplied by a factor which 
converts it to comply with the defining units used in the experiment. It 
will be observed that in all three cases the accuracy of the theory is 
adequate. Note that the theory is capable of showing the reversal in 
trend of the damping coefficient that occurs at high subsonic Mach numbers; 
however, the reversal occurs sooner and more abruptly than predicted. 

3. DAMPING IN PITCH IN THE PRESENCE OF NONLINEAR STEADY -STATE FORCES 
AND MOMENTS 

The second problem mentioned previously will now be considered; 
namely, the estimation of pitch-damping coefficients when the aerodynamic 
forces and moments are nonlinear functions of angle of attack.’ The 
procedure will again be to formulate the problem to relate as closely as 
possible to conditions in steady flow. Evidence of the existence of such 
a relationship between steady and unsteady flow can be found in the liter- 
ature; several experimental investigators (e.g., (8), ( 9 )) have noted that 
when the static stability derivative Cn^ is a nonlinear function of 
angle of attack, the pitch -damping derivative C^ + C^ is also nonlinear . 
More specifically, as indicated in Fig. 6, an increase (decrease) in 
static stability is accompanied by a decrease (increase) in pitch damping. 
Experimental results for several radically different configurations have 
confirmed the general nature of this relationship and the results, when 
cross plotted, in most cases have exhibited the linear correlation also 
shown in Fig. 6 and expressed as Eq. (l4): 

Cmq + Cny^ = A + BQn^ ( l4) 
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However, this information by itself is insufficient for the purpose of 
developing a rational estimation procedure, and the physical basis for 
such a relationship must be established. 

3.1 The Indicial Response 

In the analysis to follow, the physical basis for the relationship 
between pitch-damping and static stability will be illustrated through 
use of the concept of the indicial response. The indicial response of 
a system, by definition, is the response of the system to a disturbance 
in the form of a step function. As an illustration. Fig. 7 shows a 
step change in angle of attack Aa. The indicial response which will 
be considered, also shown, is the resulting moment coefficient variation 
divided by the increment in angle of attack. It will be noted that the 
indicial response is not zero at time zero, but has a starting value 
which may be computed rather simply (e.g., (lO)) from considerations of 
the initial momentum change imparted to the air immediately adjacent to 
the surface of the aerodynamic body. The response then proceeds to 
change with time, eventually reaching its steady-state value which has 
been indicated to be a function of angle of attack. 

The following assumptions will be made concerning the indicial 
response: First, that the total response can be divided logically into 

two individual responses as indicated in the last sketch of Fig. 7 * 

The first response is associated with the starting value, and falls to 
zero as the finite extent of the body alleviates the momentum initially 
imparted to the air. The second response is associated with the develop- 
ment of the flow pattern leading to the steady-state value. The starting 
value will be assumed independent of angle of attack. This assumption 
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is consistent vith the results of linearized (10) and second-order 
aerodynamic theory for planar bodies (wings) and will be assumed to 
hold in general. Finally, the form of the response associated wit h the 
nonlinear steady-state value will be assumed to be independent of the 
nonlinearity. This assumption is based on the intuitive reasoning that 
the propagation of flow disturbances, that is, the mechanism by means of 
which the steady-state loading is attained, is independent of the form of 
the loading. The importance of the assumption lies in the fact that it 
retains the characteristic of the response of a linear system which permits 
the principle of superposition to be applied ( ll) . With these assumptions 
the indicial response may be written as 

(t,a) = (O)fx(t) + ~ («,a)f 2 (t) (15) 

Art, Aa Ajx 

where the limits of the functions fi and f 2 are zero and unity. 

3.2 The Superposition Integral 

In Fig. 8 an arbitrary angle -of -attack variation is shown, beginning 
with an initial value a* which has existed for an infinite period prior 
to time "zero.” If the continuous variation is approximated by a series 
of step functions, then associated with each step is an indicial response 
such as previously discussed. The resulting moment variation is then 
given by the sum of the indicial responses, 

Cjit) = CmC 0 ) + (t-t 1 ,ai)[a 2 -ai] + (t-t 2 ,as) [ctQ-ae] + . . . 

= Cm(0) (t-tj^OnWn 

n 

where each response is evaluated at the appropriate instant of time (t-t n ) 
relative to its inception. Passing to the limit of an infinitesimal step 
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size then leads to Eq. ( IT ) : 

C m (t) = Qm(o) + f C™ (t-r,a(r)) dr ( 17 ) 

J o dT • 

Again one 'should note that the use of the superposition integral is possible 

only as a result of assuming an indicial response of unique form, containing 

the non lin earity associated with a only as a multiplicative factor 

determining the magnitude of the response. 

3.3 Determination of and 0 ^ 

As a result of the preceding development, it is possible to consider 
separately the two parameters and C^. which make up the total pitch 

damping. Attention will be focused here on Cn^, the pitching -moment 
coefficient proportional to vertical acceleration. The motion which leads 
to a moment proportional to a but not to q is shown in Fig. 9 - harmonic 
plunging oscillations of an aerodynamic body at constant inclination to 
the mean flight path. The angle -of -attack variation is then given by 

a(t) =051! + ctosin cat (l 8 ) 

a harmonic variation of angle of attack of amplitude <Xq about a mean 
value Ojjj. If now the steady-state value of the indicial response is 
expanded in a Taylor series about the mean angle of attack <%, 

oc) = (c%i)[ a “ c %i] + 2 cL(y2 + . • . (19) 

all the necessary equations have been developed for computing the total 
moment response to this motion. Substitution of Eqs* (l 5 ) and- (19) in 
Eq. (17) and a little manipulation yields 
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vhere 


AC m (t) = Omit) - CmCo) = 0^(0) jf^t-r) dr 

d f\ z (t-r) to - ( - T - dr + . . . 

4, dT 


+ a 


( 20 ) 


a — Qn^(®m) " 


^ (<%) + ^ (<%) + 


da 


2 da' 1 




4 - i £%* (am) ♦ • • 
6 da 


and AC^t) is the time-varying moment caused only by a time variation 
of angle of attack a. Any mean value of the moment has been subtracted 
as shown. If now the angle-of -attack variation given by Eg. (l8) is 
introduced, AC^Ct) can be determined by a straightforward integration of 
Eg. (20). It is evident that because of the nonlinearity . of the steady- 
state response, the response AC^t) to the harmonic angle of attack will 
not itself be harmonic, but will be instead a periodic function made up 
of components exhibiting some range of freguencies. However, only those 
components are of interest which are characterized by the freguency 10 , 
since only those components are capable of doing work, or of being 
interpreted as an aerodynamic spring or aerodynamic mass. Application of 
the Fourier transform to AC^t) yields the desired components. A some- 
what simpler and more direct approach is the application of the Laplace 
transform (ll) to Eg. (20); reduction to the Fourier variable, w, then 
yields the desired components directly. Either approach is straightforward 
and need not be detailed here. The final result is 
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F[ACm(t)] 
F[a(t) - a m ] 


= + iw[C I% (0)F 1 (u)) - C^FsCw)] 


( 21 ) 


where 




( 22 ) 


c mae “ 

Fi(^) = f fi(t)e iw ^dt 

J o 

/ OO 

[l-f 2 (t)]e iwt dt 

• " 

The first term in Eq. (2l) is a moment in phase with oscillatory angle of 
attack, and is therefore equivalent to an effective aerodynamic spring; 
hence the s imi larity, in the notation to the usual designation of aero- 
dynamic spring, Cj^. The second term is complex through the Fourier trans- 
forms Fi(w) and F 3 (u>). The imaginary part is in phase, with the rate , of 
change of angle of attack, and is therefore proportional to C^. The real 
part is in phase with acceleration and can thus be identified as an aero- 
dynamic mass. To the first order in frequency the aerodynamic mass 
contribution to the moment is zero. Also to the first order in 
frequency, can be expressed as 

Cn£ = C + DC^ (23) 

where 

rioo - .... 

c = ^ cw( 0) / fi(t)dt ■ • . 

1 ro ■ 

D =. - 1 F [-1 - f 2 (t)]dt 
IJ 0 

The derivative Cn^ is thus seen to have the linear correlation form 
originally given by Eq. (l4), but the correlation relates to the effective 
value rather than the static value C^. It will be seen later' 
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that this form is the correct- one. The constants C and D are functions 
of the starting value of the indicial response and integrals of the 
functions fi(t) and f 2 (t) which determine the form of the indicial 
response. The integrals are equal to the areas shown in Fig. 10. 

The reader will note that the effective static stability Cmae gl ven 
by Eq. (22) is the value which would be obtained by means of an experimental 
apparatus wherein the motion is inexorably forced and the aerodynamic 
moment component in phase with the angle of attack is measured. It can 
also be shown by the Kryloff-Bogoliuboff method (12) of nonlinear mechanics 
to be the value which would be obtained from frequency measurements of a 
tuned apparatus, that is, an apparatus which is forced to oscillate at 
its resonant frequency. It was by the latter method that the data to be 
presented subsequently were obtained. 

A similar analysis for the derivative Cm^ leads to the result 

Qmq = Cmq^(ct(m) (2^-) 

which states that, to first order in frequency, Cm^ is equal to the 
effective steady-state value Cm^ existing at the mean angle of attack <%. 

3.k Discussion 

From the previous development it is evident that the total pitch 
damping, as given by the sum of Eqs. ( 23 ) and (24), 

CjDq + = Cmq^ + C + DCm^ (25) 

should exhibit a linear correlation only if Cm^ is independent of angle 
of attack, varies with angle of attack in the same manner as Cm^, or is 
negligibly small. It is noted that those configurations investigated to 
date which have failed to exhibit a linear correlation have all been winged 
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configurations which would be expected to have a sizable value of 
relative to . C^. Bodies of revolution, which would be expected to have 
a much smaller value of C^, have all correlated rather well. A typical 
example is shown in Fig. 11 where the pitch danqping and the effective value 
of static stability are shown as functions of mean angle of attack for 
three oscillation amplitudes. These are data recently obtained at Ames 
Research Center for a body of revolution at a subsonic Mach number. The 
value of Cm^ for this configuration is estimated to be quite small and, 
as can be seen, the data correlate very well. Rote that, particularly in 
the region 'of greatest nonlinearity, the pitch damping varies with oscil- 
lation amplitude and would not, therefore, correlate with the single static 
stability curve given by C^. It is the simultaneous variation of the 
effective static stability Ci^ which leads to the correlation shown. 

It appears, therefore, that a fairly sound understanding of the 
physical processes underlying a correlation between the dynamic stability 
derivative C^ and the static stability derivative C^ has been 
attained. The constants in Eq. (23). indicate the quantities which must be 
determined in order to utilize these results in a rational estimation 
procedure - the starting value of the indicial response, and the charac- 
teristics of the response as determined by integrals of the functions fi(t) 
and f 2 (t). Work in this direction is now in progress. 
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Fig. 1.- Variation of 
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®L and x^, vith reduced aspect ratio, PA. 
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Fig. 2.- Variation of lift and moment coefficient corection terms vith 

reduced aspect ratio, PA. 




Fig. 3.- Comparison of theoretical and experimental damping coefficients 
for a triangular wing of aspect ratio 1.45* 
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Fig. 4.- Comparison of theoretical and experimental damping coefficients 
for a triangular wing of aspect ratio 2. 
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• - Comparison of theoretical and experimental damping coefficients 
for a triangular ving of aspect ratio 4. 
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Fig. 6.- Relationship between static stability and damping in pitch. 
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Fig. 7*~ The indicial response 



Fig. 8.- Approximation of continuous angle -of -attack variation by means 

of step functions . 


DIRECTION OF MOTION 



Fig. 9 .- Motion which produces harmonic angle -of -attack variation. 



Fig. 10.- Interpretation of the integrals appearing in the correlation 

constants. 



Fig. 11.- Experimental correlation of effective static stability and 

damping in pitch. 


